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We perform radiation-hydrodynamics simulations of binary neutron star mergers in numerical 
relativity on the Japanese “K” supercomputer, taking into account neutrino cooling and heating 
by an updated leakage-plus-transfer scheme for the first time. Neutron stars are modeled by three 
modern finite-temperature equations of state (EOS) developed by Hempel and his collaborators. We 
find that the properties of the dynamical ejecta of the merger such as total mass, average electron 
fraction, and thermal energy depend strongly on the EOS. Only for a soft EOS (the so-called SFHo), 
the ejecta mass exceeds O.OIMq. In this case, the distribution of the electron fraction of the ejecta 
becomes broad due to the shock heating during the merger. These properties are well-suited for 
the production of the solar-like r-process abundance. For the other stiff EOS (DD2 and TMl), for 
which a long-lived massive neutron star is formed after the merger, the ejecta mass is smaller than 
O.OIM 0 , although broad electron-fraction distributions are achieved by the positron capture and 
the neutrino heating. 

PACS numbers: 04.25.D-, 04.30.-w, 04.40.Dg 


I. INTRODUCTION 

The merger of binary neutron stars (BNS) is one of 
the most promising sources of gravitational waves for ad¬ 
vanced LIGO [I| , advanced VIRGO 0 , and KAGRA Q , 
which will start operation in a few years. The recent 
statistical studies suggest that these gravitational-wave 
detectors will observe gravitational waves from merger 
events as frequently as ~ 1-100/yr ii- The merger of 
BNS is also a promising candidate for the central engine 
of short-hard gamma-ray bursts. If gravitational waves 
are observed simultaneously with them, a long-standing 
puzzle on the central engine of short-hard gamma-ray 
bursts may be resolved. 

In addition to these aspects, BNS are attracting at¬ 
tentions as the nucleosynthesis site of heavy elements 
by the r-process Q, which may proceed in the neutron- 
rich matter ejected during the merger. Recent obser¬ 
vations of metal-poor stars Q strongly suggest that 
there should exist ’main’ r-stars affected by ’universal’ r- 
process events in which the resulting abundance is close 
to that of solar-abundance pattern for nuclei with the 
atomic number ^ ^ 38 (A > 90). It has recently been 
revealed ii that the supernova explosion, which was 
previously considered to be the most promising candidate 
for the site of the r-process, may not be a viable origin 
in this regard, and the BNS mergers is getting attention. 

Furthermore, a strong electromagnetic emission may 
accompany the radioactive decay of the r-process ele¬ 
ments dMl and it could be an electromagnetic coun¬ 
terpart of gravitational waves from BNS mergers. An 
infrared transient event associated with GRB 130603B 
is the first candidate for such events [l^. These facts 
strongly encourage the community of gravitational-wave 
astronomy to explore the r-process nucleosynthesis and 


associated electromagnetic emission in the BNS merger. 

For the quantitative study of these topics, we have to 
clarify the merger dynamics, subsequent mass ejection, 
and physical condition of the ejecta, which are neces¬ 
sary to study the nucleosynthesis, subsequent decay of 
the heavy elements in the ejecta, and electromagnetic 
emission from the ejecta. For this purpose, we have 
to perform BNS merger simulations taking into account 
both general relativistic gravity and detailed microphys¬ 
ical processes. 

For the former, recent numerical relativity simulations 
(e.g., [l^; see also [l^ for simulations in approximate 
general relativistic gravity) have clarified that the gen¬ 
eral relativistic gravity can be the key for the mass ejec¬ 
tion: In general relativity, shock heating plays a promi¬ 
nent role in the merger process, and consequently, the 
ejecta that is dynamically expelled during the merger 
(dynamical ejecta) are composed not only of those driven 
by the tidal interactions but of those driven by the ther¬ 
mal pressure, by contrast with the result in Newtonian 
simulations (e.g., El) for which the tidal component is 
major. 

For the latter, we have recently developed a neutrino- 
radiation hydrodynamics code, and now, we can perform 
simulations both employing a wide variety of equations 
of state (EOS) for the nuclear matter in which finite- 
temperature effects are incorporated and handling neu¬ 
trino cooling and heating with reasonable sophistication. 

This is the first study based on these modern aspects 
of the merger dynamics in general relativity taking into 
account the microphysics. In this paper, we report the 
latest result of our simulations for equal-mass BNS merg¬ 
ers of typical neutron-star mass (I. 35 M 0 ) for three repre¬ 
sentative EOS, among which the radius of neutron stars 
is appreciably different. In this paper, we only consider 
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EOS no Eg (MeV) K (MeV) g (MeV) L (MeV) 


SFHo 

0.1583 

16.19 

245.5 

31.57 

47.10 

DD2 

0.1491 

16.02 

242.7 

31.67 

55.03 

TMl 

0.145 

16.3 

281 

36.9 

110.8 


TABLE I. Characteristic properties of EOS at the nuclear 
saturation density, no: the nuclear saturation density. Eq: 
the binding energy. K : the incompressibility. S : the sym¬ 
metry energy. L : the logarithmic derivative of the symmetry 
energy. 

the case of equal-mass binaries. The dependence on the 
mass-ratio and the total mass will be studied in a future 
work. We will show that the physical properties of the 
dynamical ejecta such as the mass and neutron fraction 
depend strongly on the EOS. We find that for producing 
mildly neutron-rich dynamical ejecta of large mass with a 
broad range of the electron fraction, a relatively soft EOS 
that yields small-radius (< 12 km) neutron stars is neces¬ 
sary. Because of such a broad distribution of the electron 
fraction, the universal Q solar-abundance pattern of the 
r-process elements may be reproduced without need for 
the other contributions M- 


II. METHOD, EOS, INITIAL MODELS, AND 
GRID SETUP 

We solve Einstein’s equation by the puncture-BSSN 
(Baum gar te-Shapiro-Shibata-Nakamura) formalism as 
before [ll, [l^ . The 4th-order finite-differencing scheme 
is applied to discretize the field equations. The radia¬ 
tion hydrodynamics equations are solved by a recently- 
developed code which is updated from the previous ver¬ 
sion: In this new code, neutrino transport is computed 
in a leakage-based scheme incorporating Thorne’s 
moment formalism with a closure relation for a free- 
streaming component [2lj. For neutrino heating, absorp¬ 
tion on free nucleons is taken into account. 

We employ three EOS for nuclear matter derived re¬ 
cently by Hempel and his collaborators, which are re¬ 
ferred to as SFHo [13], DD2 [H, and TMl [H in the fol¬ 
lowing. TMl EOS, which is also known as Shen EOS [^ . 
is based on the relativistic mean field theory with a pa¬ 
rameter set of Ref. [2^ and have been used widely in 
both supernova and compact-binary merger simulations. 
SFHo EOS is constructed so that the predicted neutron 
star radius matches recent neutron star observations by 
extending the nonlinear Walecka model [131 . DD2 EOS 
is based on a relativistic mean field model with a density 
dependent coupling [131 . Some characteristic properties 
of EOS are listed in Table H) 

For all of them, the predicted maximum mass for 
spherical neutron stars is larger than the largest well- 
measured mass of neutron stars, ss 2Mq [1^. For these 
EOS, the radius of neutron stars with mass 1.35M0 
is Ri .35 = 11.9 km (SFHo), 13.2 km (DD2), 14.5 km 


Model 

Ri.3S (km) 

Axg (m) 

N 

Mej (Mo) 

(Ue) 

SFHo (high) 

11.9 

150 

285 

1.1 X lO”’^ 

0.31 

SFHo (low) 


250 

160 

1.3 X 10"^ 

0.32 

SFHo (no-heat) 


250 

160 

1.0 X lO”’^ 

0.29 

DD2 (high) 

13.2 

160 

285 

2.1 X 10-3 

0.29 

DD2 (low) 


270 

160 

1.9 X 10“3 

0.29 

DD2 (no-heat) 


270 

160 

0.9 X 10-3 

0.26 

TMl (high) 

14.5 

200 

285 

1.2 X 10-3 

0.26 

TMl (low) 


300 

160 

0.8 X 10-3 

0.25 


TABLE 11. i?i, 35 : the radius of spherical neutron stars of 
mass 1.35M0. Axg: the grid spacing in the finest refinement 
level. N: the grid number in one positive direction for each 
refinement level. Mej and {Ys) denote the ejecta mass and the 
averaged value of Ye measured at the end of the simulations. 
Model name follows the EOS. 

(TMl), respectively (see Table HIl) . We refer to an EOS 
with a small neutron star radius {R 1.35 < 12 km) like 
SFHo as a soft EOS and an EOS with a large radius 
{R 1.35 > 13 km) as a stiff EOS. The stellar radius plays 
a key role for determining the merger remnant and the 
properties of the dynamical ejecta. 

In numerical simulations, we have to follow the ejecta 
with velocity 0.1-0.3c (c is the speed of light), which ex¬ 
pand to > 10^ km in the simulation time. To follow the 
ejecta motion as well as to resolve neutron stars, we em¬ 
ploy a fixed mesh-refinement algorithm. In this work, we 
prepare 9 refinement levels with the varying grid spacing 
as Axi = 2®“* Axg {I = 1,2, ■ ■ ■ ,9) and all the refinement 
levels have the same coordinate origin. Here, Axi is the 
grid spacing for the l-th level in the Cartesian coordi¬ 
nates. For each level, the computational domain covers 
the region [—NAxi,NAxi] for x- and y-directions, and 
[0, AAa;;] for z-direction (the reflection symmetry with 
respect to z = 0 is imposed). In the highest-resolution 
run, we assign N = 285, Axg = 150-200 m, and utilize 
« 7, 000 CPUs on the K computer. 

To check that the numerical results depend only weakly 
on the grid resolution, we also performed lower-resolution 
simulations. For this case, N = 160 and Aig = 250- 
300 m. As listed in Table im we found that the results 
such as total ejecta mass and averaged values of Ye de¬ 
pend very weakly on the grid resolution. Furthermore, to 
confirm the importance of the neutrino heating, we also 
performed simulations in which the neutrino absorption 
is switched off (denoted as ’no-heat’ in Table HIl) and com¬ 
pared the results for the first time. 

We consider equal-mass BNS with each mass 1.35M0. 
Observed neutron stars in BNS typically have the mass 
ratio close to unity and the mass in the range 1.20- 
I. 45 M 0 [33]. Thus, our choice reasonably reflects the 
observational fact. The initial orbital separation is cho¬ 
sen so that the orbital angular velocity, U, satisfies 
GmoEl/c^ = 0.028 where mg = 2.7Mq is the sum of each 
mass in isolation and G gravitational constant, respec¬ 
tively. Table m lists the key parameters of our models 
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FIG. 1. Mass (upper panel) and characteristic velocity (lower 
panel) of the ejecta as functions of time for SFHo (red solid), 
DD2 (blue dashed), and TMl (green dotted-dashed), tjvf-e 
approximately denotes the time at the onset of merger (see 
text). 


and simulation setup. 


III. RESULTS 

For all the models, a massive neutron star (MNS) is 
formed after the onset of merger as expected from our 
previous results [s^. The MNS are long-lived in the 
sense that their lifetime is much longer than their ro¬ 
tation period of < 1 ms. For SFHo, the MNS eventually 
collapses to a black hole (BH) in ~ 10 ms because the 
maximum mass of spherical neutron stars is relatively 
small as ~ 2.0Mq. The mass and spin parameter of the 
BH are Mbh ~ 2.6Mq and obh ~ 0.70, and a torus 
with mass Mtorus ~ 0.05Mq is formed around it. Such 
a system may be a central engine of short-hard gamma- 
ray bursts. For other two cases, the remnant MNS does 
not collapse to a BH in our simulation time ~ 30-40 ms. 
Because the maximum mass of spherical neutron stars 
for DD2 and TMl is « 2.4 and 2.2Mq, the formed hot 
and rapidly rotating MNS with mass 2.6Mq will not 
collapse to a BH unless a substantial fraction of the angu¬ 
lar momentum and thermal energy is dissipated by some 
transport process and the neutrino emission, respectively 
(e.g., [il 0 ). 

Figure [T] plots the evolution of the rest mass Mej and 
the characteristic velocity Kj for the ejecta. Here, Im-g 
denotes the time at which Mej exceeds 10~^Mq (here¬ 
after we will use tM-e as the time at the onset of merger). 
We specify the matter as the ejecta if the time component 
of the fluid four velocity ut is smaller than —1. Note that 
another condition [3l| for the ejecta hut < — 1 where h 
is the specific enthalpy, which may be more appropriate 


for the hot matter, gives slightly larger ejecta mass. 14 j 
is defined by y^ 2 £'kin/Mej where ifkin is kinetic energy of 
the ejecta. Figure [T] shows that the ejecta mass depends 
strongly on the EOS: For softer EOS (i.e., for smaller 
values of i^i.as), the ejecta mass is larger. Remarkable 
is that with the decrease of i?i ,35 by ~ 3 km, the ejecta 
mass increases by more than one order of magnitude and 
only for R 1.35 < 12km the ejecta mass exceeds O.OIMq, 
as already indicated in [l3^[l3■ The av erag ed ejecta ve¬ 
locity is ~ 0.1-0.2c as also found in (l^. [l5l|. In the later 
phase, the total ejecta mass relaxes approximately to a 
constant, and the ejecta are in a free expansion phase for 
all the models. 

There are two major mass ejection mechanisms during 
the merger phase. One is tidal interaction and the other 
is shock heating. By the tidal interaction, the matter 
tends to be ejected near the orbital plane. On the other 
hand, by the shock heating, the matter is ejected in a 
quasi-spherical manner. Because both effects play a role, 
the ejecta usually have a spheroidal morphology. For 
small values of R 1 . 35 , the shock heating plays a stronger 
role and the ejecta in this case have a quasi-spherical 
morphology. 

Figured] plots the profiles of the electron fraction, Y^, 
(left half) and entropy per baryon, s, (right half) of the 
ejecta on the x-y and x-z planes for DD2 (left panel) and 
SFHo (middle and right panels). For DD2, the ejecta are 
composed of (i) tidally-ejected matter with low values of 
Yg and s near the orbital plane and (ii) shock-heated mat¬ 
ter with relatively high values of Yg. The shock-heated 
ejecta are less neutron-rich because the temperature gets 
much higher than ^ 1 MeV as a result of the shock heat¬ 
ing, producing copious pairs that activate e~ and 

e"*" captures by protons and neutrons, respectively. As a 
result of e~ and e"*" captures, the luminosities of Vg and Dg 
become quite high as > 10®^ ergs/s (see Fig. [3]), as long as 
the remnant MNS is present. Because the original ejecta 
are neutron-rich, e’*' capture dominates e~ capture, and 
hence, the luminosity of Dg is higher than that of iZg [T^ 
and the ejecta become less neutron-rich. 

In addition to the tidal-driven and shock-heated com¬ 
ponents explained above, we found the third compo¬ 
nent in a later phase, that is, neutrino-heated compo¬ 
nent with even higher values of Yg and s in the re¬ 
gion above the MNS pole (see the high-entropy region 
in the left panel {x-z plot) of Fig. [3]). Furthermore, 
some fraction of the material obtains enough energy to 
be additional neutrino-driven ejecta. Possible existence 
of such a component was recently reported in a MNS sys¬ 
tem and a BH and torus system which is expected 

to be formed after the BNS mergers [s^. We confirmed 
the existence of the neutrino-driven component in self- 
consistent numerical-relativity simulations of the merger 
for the first time. 

For TMl, the results are basically similar to those for 
DD2 except for the fact that the tidally-ejected compo¬ 
nent is more dominant and the e'^ capture is less effi¬ 
cient. Also, the neutrino-driven wind appears to play 
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FIG. 2. Contours of the electron fraction, Ye, (left half) and the entropy per baryon, s, (right half) in x-y (lower) and x-z 
(upper) planes, left panel: for DD2 at 8.5 ms after the merger, middle panel: for SFHo at 5.0 ms after the merger, right panel: 
for SFHo at 15.0 ms after the merger. 


a major role for the mass ejection (see the curve for 
t — tM-6 > 5 ms of Fig. [U because the total ejecta mass 
for this EOS is rather small. Here, note that it is not 
easy to exclude the effect of artificial atmosphere in grid- 
based simulations, in particular when the ejecta mass is 
low (< 10 “^Mq) as in the case of TMl. The contam¬ 
ination in mass would be ^ when the ejecta 

expand to ^ 2000 km in our setting of the atmosphere 
with density ^ lO^g/cm^, while it would be of order of 
percent if the ejecta is as massive as ^ The 

contamination in Yg would be similar level. For this rea¬ 
son, in the following, we will basically consider DD2 as a 
representative of a stiff (or moderately stiff) EOS. 

For SFHo, shock waves are formed for several times 
during the merger phase as the MNS oscillates with a 
high amplitude, and hence, a certain fraction of mat¬ 
ter originally ejected by the tidal interaction is subse¬ 
quently heated up by shocks (s increases), resulting in 
the increase of the values of Ye via weak interactions. On 
the other hand, other parts less influenced by the shock 
heating preserve the neutron-rich nature of the original 
neutron stars. As a result of these two facts, the ejecta 
can have higher values of s and Ye than for DD2 and 
TMl even in the orbital plane with an appreciably in¬ 
homogeneous distribution of Ye (see the middle panel of 
Fig. [2|). Because a BH is formed at ^ 10 ms after the 
onset of merger for SFHo, the strong neutrino emission 
region is swallowed into the BH and neutrino luminos¬ 
ity decreases to < 10®^ ergs/s. Hence, there is less clear 


neutrino-driven ejecta component for this EOS (see the 
bottom panel of Fig. Ej). 

The upper panel of Fig. E] shows the time evolution of 
averaged values of Ye ((Fe)) from which the effect on Ye of 
the shock heating and the resulting positron capture can 
be seen more clearly. The several distinct changes in (Yg) 
observed for SFHo in < Sms after the onset of merger 
reflect the strong e+ capture activated by the shock heat¬ 
ing. During this phase, (Ye) for SFHo increases drasti¬ 
cally to be « 0.3. After this phase, on the other hand, 
(Ye) for SFHo is approximately constant because the e~ 
and e"*" captures balances and because the neutrino lu¬ 
minosity decreases to be ~ 10®^ ergs/s due to the BH 
formation, which is not sufficient to change (Ye) of the 
massive ejecta. Thus, for softer EOS like SFHo, Yg is 
likely to be increased primarily by the e+ capture. 

On the other hand, (Yg) for DD2 and TMl in the early 
stage is low as Yg < 0.1-0.2, while it increases in time. 
This is simply because the shock heating at the first con¬ 
tact is not strong enough to increase (Yg) significantly 
for these stiffer EOS; i.e, the original composition of the 
ejecta driven by tidal torque, which is composed primar¬ 
ily of neutron-rich matter with low temperature, is tem¬ 
porally preserved as found in 0,[3. In the later phase, 
however, the ejecta become less neutron-rich. This is 
partly due to the positron capture discussed above. In 
addition, the electron neutrinos emitted from the rem¬ 
nant MNS convert some fraction of neutrons to protons 
via the electron neutrino capture (see below for a more 
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FIG. 3. Luminosity curves of (red solid), zJe (blue dashed), 
and heavy (green dotted-dashed) neutrinos for TMl (top), 
DD2 (middle), and SFHo (bottom). 


detailed discussion). For stiffer EOS, the importance 
of the electron neutrino capture in increasing Yf. of the 
ejecta is enhanced because of their lower temperature 
and the maintained high neutrino luminosity from the 
long-lived MNS. 

The lower panel of Fig. 2] plots the mass-distribution 
histograms for Fg normalized by the total mass of the 
ejecta at « 25 ms after the onset of merger. For all of 
the models, is distributed in a broad range between 

0.05 and 0.45. This result is com plet ely different from 
that found in the previous studies [1^ [l^ in which the 
distribution of Y^ is very narrow with a lower average 
value <0.1. This disparity can be explained as follows. 

In the previous approximate general relativistic 
study [l^ , the weak interaction processes were not taken 
into account, and hence, the ejecta remain neutron-rich 
because there is no way to change Y^. In the previous 
Newtonian studies [l^, they took into account the neu¬ 
trino cooling (e“ and captures). However, as we men¬ 
tioned already, the effect of the shock heating is under¬ 
estimated significantly in Newtonian gravity, and hence, 
the effect of the capture would be much weaker than 
that in our simulations due to the underestimated tem¬ 
perature. In addition, they did not take into account the 
neutrino heating (absorptions) which is expected to play 
a role for stiffer EOS in which the positron capture is 
relatively less important due to lower temperature. 

To see the effects of the neutrino heating more quan¬ 
titatively, we performed simulations without (no-heat) 
neutrino heating for SFHo and DD2. We found that for 
both EOS, the contribution of the neutrino-driven com¬ 
ponent in the ejecta mass is ^ 1Q~^Mq at the end of the 
simulation (see Table El), which is consistent with that 


FIG. 4. Upper panel: The time evolution of the averaged 
value of Ye for SFHo (red solid), DD2 (blue dashed), and TMl 
(green dotted-dashed). Lower panel: The mass-distribution 
histograms of Ye normalized by the total mass of ejecta mea¬ 
sured at ~ 25 ms after the onset of merger for SFHo, DD2, 
and TMl. 

found in The amount of the neutrino-driven ejecta 
is minor for SFHo but comparable to the amount of the 
dynamical ejecta for DD2. This result suggests that the 
neutrino heating plays a relatively more important role 
for stiffer EOS like DD2 and TMl in which the amount 
of the dynamical ejecta is 10 “^Mq. 

The neutrino heating plays an important role in chang¬ 
ing the chemical composition (Ye) of the ejecta. As shown 
in Fig. [31 the luminosities of Vg and Vg are quite high as 
> 10®^ ergs/s. Due to the absorption of neutrinos with 
this high luminosity, the ejecta become more proton-rich 
because the electron neutrinos convert some fraction of 
neutrons to protons via the reactions n + Vg p + e~. 
Note again that Vg capture is more efficient than Vg cap¬ 
ture since the ejecta are neutron-rich. 

Figure [S] compares the time evolution of (Yg) (upper 
panel) and the mass-distribution histograms for Yg at 
Ri 25 ms after the onset of merger (lower panel) between 
simulations with and without neutrino heating for SFHo 
and DD2. The results indicate that for SFHo, (Yg) is in¬ 
creased to be Ri 0.29 due to the positron capture and the 
neutrino heating pushes up it further by Ri 0.02 at the end 
of the simulations. For DD2, the effect of the positron 
capture is weaker and the neutrino heating plays a rel¬ 
atively important role, increasing (Yg) by r; 0.03. Such 
enhancements of (Yg) due to the neutrino heating would 
be important in considering the r-process nucleosynthe¬ 
sis [13. 

The mass-distribution histograms also shift towards 
the higher Yg side due to the neutrino heating. How¬ 
ever, the distributions still show a broad feature even 
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Ye 

FIG. 5. Same as Fig.[4]but for simulations with and without 
(denoted as no-heat) the neutrino heating for SFHo (red and 
magenta (no-heat)) and DD2 (blue and light blue (no-heat)). 


without the neutrino heating. This suggests that the 
positron capture resulting from the strong shock heating 
due to general relativistic gravity is primarily responsi¬ 
ble for making the Ye distribution broad for DD2 and 
SFHo. For much stiffer EOS like TMl, the neutrino 
heating would play a relatively major role. Although 
our treatment for the neutrino transfer is an approximate 
one, our results indicate that the neutrino heating plays 
an important role in determining the chemical properties 
of the ejecta. 


IV. SUMMARY AND DISCUSSION 

We have reported the hrst numerical results of radia¬ 
tion hydrodynamics simulations in general relativity fo¬ 
cusing on the properties of the dynamical ejecta of the 
equal-mass BNS merger with typical mass of each neu¬ 
tron star (1.35Mq). Three modern finite-temperature 
EOS are employed to clarify the dependence of the ejecta 
properties on the EOS. We found that the total mass 
of the ejecta is larger for softer EOS (giving smaller- 
radius neutron stars), and it exceeds O.OIMq only for 
the case that R 1.35 < 12 km, as indicated in H- As 
shown in [13, [13 1 the electromagnetic luminosity of the 
ejecta by the radioactive decay of the r-process elements 
would depend sensitively on the ejecta mass, and hence, 
the predicted range of the luminosity spans in a wide 
range due to the uncertainty of the nuclear-matter EOS. 

We also found that the averaged value of Ye of the 
ejecta is higher for softer EOS like SEHo in which i?i ,35 


is smaller, reflecting the fact that the shock heating is 
more efficient. For all of the models, the value of Ye for 
the ejecta has a broad distribution between ^ 0.1 and 
0.45, by contrast with the previous studies [il,[i3. Here, 
both the strong shock associated with general relativistic 
gravity and the weak interactions play crucial roles for 
this. Such a broad distribution may be well-suited for 
producing the universal Q solar-abundance pattern of 
r-process elements as illustrated in [i3- 

For the EOS but for SFHo, the dynamical ejecta mass 
is of order 10 “^Mq. In this case, a rather higher merger 
rate of > 10“"^ yr~^ than the present estimates of the 
Galactic rate (a few 10“^ yr~^) HI] is necessary to ex¬ 
plain the amount of heavy r-process elements [Sa, 1^ , if 
the the dynamical ejecta from binary neutron star merg¬ 
ers is responsible for their production. In regards to this 
point, SFHo is an attractive EOS. We will study conse¬ 
quences of our results on the synthesis of heavy elements 
in the forthcoming paper. If EOS is not very soft like 
SEHo, some other contributions, such as mergers of black 
hole-neutron star binaries , disk winds from accretion 
torus around a merger remnant black hole [13,113 , and 
magnetorotational supernova explosions [doj ] may be nec¬ 
essary. In such cases, however, it is not clear whether the 
universality requirement can be achieved or not. 

In this work, we focused only on the equal-mass bi¬ 
nary case and did not explore the dependence of the re¬ 
sults on the binary parameters such as the total mass 
and the mass ratio. As reported in [l3| . the relative im¬ 
portance of the tidal interactions and the shock heating 
in the dynamical mass ejection depends on the binary 
parameters. It is interesting to explore the dependence 
of the results on binary parameters for SEHo and the 
resulting abundance profile in the future work, because 
the observed abundance patterns of the metal-poor, r- 
rich stars show some diversity in the lower mass-number 
region Q. Also, we did not continue our simulations be¬ 
yond 30-40 ms after the onset of merger. Eor the longer 
time scales, magnetohydrodynamic processes (4lj |. vis¬ 
cous heating, and nuclear recombination [13 could be 
important. Self-consistent studies of these effects in the 
BNS merger also have to be done in the future. 
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